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1.0  INTRODUCTION 


From  March  1980  to  March  1983  the  University  of 
Lowell  Center  for  Atmospheric  Research  (ULCAR)  conducted  iono¬ 
spheric  radio  research  using  Digisonde  data  in  support  of 
AFGL's  ongoing  research  work  in  the  polar  cap,  the  auroral  re¬ 
gions  and  at  the  magnetic  equator.  The  complicated  structures 
of  the  ionization  in  these  regions  result  in  different  signa¬ 
tures  of  the  HF  radio  echoes.  Interpretation  of  the  observa¬ 
tions  is  therefore  not  easy,  but  the  multi-parameter-1  frgisonde 
ionograms  helped  us  in  understanding  some  of  the  processes 
that  take  place.  This  final  report  summarizes  the  results  of 
our  observations  and  activities,  and  frequent  reference  is 
made  to  scientific  reports  and  papers  that  have  been  published 
in  the  course  of  our  work. 

Section  2  discusses  vertical  and  backscatter  iono- 
gram  observations  and  relates  them  to  aircraft  measurements 
and  to  local  magnetic  activity.  A  summary  of  the  scaled 
hourly  ionogram  values  for  Goose  Bay  for  the  period  January 
1979  to  September  1982  is  given  in  the  appendix.  Special  ob¬ 
servational  programs  in  support  of  the  over-the-horizon  radar 
tests  are  described  in  Section  3.  Section  4  summarizes  the 
results  of  the  Doppler-Dr if t  measurements  with  the  Digisonde 
128PS  at  Goose  Bay.  The  drift  data  that  were  analyzed  seem 
to  indicate  that  the  Goose  Bay  ionosphere  is  connected,  at 
least  part  of  the  time,  to  the  two-cell  polar  cap  convection 
pattern.  The  direction  of  the  drift  velocity  vectors  is  con¬ 
sistently  westward  in  the  evening  hours  and  switches  to  east¬ 
ward  at  local  midnight.  Although  not  enough  data  were  ana¬ 
lyzed  to  generalize  the  results  we  believe  now  that  HF  Doppler 
drift  measurements  are  a  powerful  and  low  cost  method  to  study 
motions  in  the  ionosphere.  Digisonde  and  optical  measurements 
in  the  polar  cap  at  Thule,  Greenland,  show  during  quiet  condi¬ 
tions  the  dawn-to-dusk  drift  of  F-region  arcs,  and  F-region 


ionization  patches  moving  anti-sunward  during  disturbed  mag¬ 
netic  conditions  (Section  5). 

One  of  the  annoying  facts  in  ionogram  based  research 
is  the  normally  tedious  process  of  reducing  the  ionograms  to 
the  vertical  electron  density  distribution.  After  years  of 
study  we  succeeded  to  develop  fully  automatic  methods  of  scal¬ 
ing  the  Digisonde  ionograms  and  calculating  the  electron  pro¬ 
files  (Section  6).  The  Automatic  Real  Time  Ionogram  Scaler 
with  True  Height  Analysis  (ARTIST),  now  operating  at  Goose  Bay 
with  the  Digisonde  128PS  is  the  first  and  only  instrument,  as 
far  as  we  know,  that  outputs  the  standard  ionospheric  param¬ 
eters  and  profiles  in  real  time,  even  under  disturbed  condi¬ 
tions.  This  breakthrough  will  make  it  possible  in  the  future 
to  instantaneously  update  the  global  ionospheric  model,  thus 
greatly  improving  short  term  HF  and  VHF  predictions. 

In  1982,  AFGL  organized  an  equatorial  spread  F  (ESF) 
campaign  including  chemical  releases  from  rockets  launched  at 
Natal,  Brazil.  The  Brazilian  Ionospheric  Modification  Exper¬ 
iments  (BIME)  were  designed  to  test  the  understanding  of  the 
nighttime  ESF  phenomenon.  ULCAR  participated  with  Digisonde 
observations  at  West  Catre,  near  Natal,  and  airborne  measure¬ 
ments.  The  Digisonde  ionograms  monitored  F-layer  height  and 
upward  velocity,  specifying  the  right  launch  time.  Incidence 
angle  and  Doppler  informations  in  the  ionograms  allowed  to 
track  the  chemical  releases  and  natural  irregularities  (Sec¬ 
tion  7). 

Airborne  Digisonde  observations  were  also  conducted 
during  another  active  ionospheric  experiment.  In  September 
1981,  the  HF  heating  facility  near  Arecibo  was  illuminating 
the  F-region  at  7  and  5  MHz.  The  airborne  sounder  recorded 
oblique  echo  returns  from  the  heated  volume  only  at  specific 
aircraft  positions.  It  is  believed  that  echoes  were  received 
when  ray  paths  at  right  angles  to  the  field  aligned  irregular¬ 
ities  existed  (Section  8). 
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2.0  IONOGRAM  OBSERVATIONS  AT  GOOSE  BAY 


The  Goose  Bay  Ionospheric  Observatory,  operated  by 
personnel  of  Canadian  Marconi  for  AFGL,  uses  a  Digisonde  128PS 
system  (Bibl  and  Reinisch,  1978)  to  record  vertical  and  back- 
scatter  ionograms. 


2.1  Observational  Programs 

At  the  present  time  four  vertical  and  four  back- 
scatter  ionograms  per  hour  are  recorded  on  magnetic  tape  (7 
track,  556  cpi).  Every  other  ionogram  pair  is  printed  on  an 
electrostatic  plotter  showing  side-by-side  the  amplitude  ion¬ 
ogram  with  numbers  in  4  dB  steps  and  the  status  ionogram  with 
numbers  expressing  Doppler  and/or  directional  information. 

The  vertical  ionograms  use  a  vertical  rhombic  antenna  for 
transmission  and  an  array  of  four  crossed-loop  antennas  for 
reception.  The  backscatter  ionograms  use  a  log-periodic 
transmitting  antenna  and  a  linear  twelve  antenna  (log-periodic 
loop)  array  for  reception.  Until  January  1983  the  boresight 
of  the  linear  array  was  magnetic  north  (32°  west  of  geographic 
north).  The  auroral  oval  is  not  centered  at  the  magnetic  pole 
and  as  a  result  the  oval  aligned  ionospheric  features  such  as 
the  trough  and  the  poleward  trough  wall  approach  in  the  even¬ 
ing  sector  from  a  direction  between  30°  east  of  magnetic  north 
(early  evening)  to  20°  east  of  magnetic  north  (later  evening). 
To  improve  the  signal-to-noise  ratio  for  backscatter  echoes 
from  the  oval  aligned  features,  we  turned  the  boresight  of 
the  array  by  25°  to  the  east,  so  that  since  mid-January  1983 
the  boresight  is  7°  west  of  geographic  north.  As  of  now  we 
have  not  verified  whether  the  intensity  and/or  frequency  of 
occurrence  of  the  backscatter  echoes  has  increased  as  we  ex¬ 
pected.  It  should  be  mentioned  that  the  log-periodic  trans¬ 
mitter  antenna  is  still  pointing  due  magnetic  north.  For  the 


vertical  ionograms  the  receiving  array  consists  of  four 
crossed-loop  antennas  that  allow  beam  steering  and  polariza¬ 
tion  measurements.  All  vertical  ionograms  compare  the  verti¬ 
cal  0  and  X  echoes  and  signals  from  two  slightly  oblique  di¬ 
rections.  It  would  require  too  much  sounding  time  to  probe 
many  different  directions  for  echoes.  The  Digisonde  allows 
azimuthal  scans  in  30°  increments  for  two  fixed  zenith  angles: 
18°  and  38°.  Until  mid-January  1983  the  hourly  ionograms 
steered  the  receiving  beam  toward  magnetic  north  at  the  two 
zenith  angles  of  18°  and  38°,  while  the  other  vertical  iono¬ 
grams  probed  the  azimuth  angles  of  ±30°  around  magnetic  north 
at  a  zenith  angle  of  approximately  18°.  No  clear  indication 
of  ionization  moving  in  the  magnetic  north-south  or  east-west 
direction  was  obtained.  So  in  January  1983,  the  beam  direc¬ 
tions  were  changed.  The  hourly  ionogram,  i.e.  the  "C  program" 
with  Z  =  1,  steers  the  beam  toward  magnetic  north  and  60°  east 
of  north  at  a  zenith  angle  of  18°  (in  addition  to  the  verti¬ 
cal  direction).  The  Z-parameter  in  the  Digisonde  128PS  se¬ 
lects  the  desired  antenna  steering  directions  (see  Table  1). 
The  other  vertical  ionograms,  i.e.  the  "A  program"  with  Z  =  9, 
looks  30°  and  210°  east  of  magnetic  north  (zenith  angle  = 

18°),  essentially  sampling  the  geographic  north-south  direc¬ 
tion.  We  have  not  yet  analyzed  the  recent  data  in  terms  of 
north-south  propagating  disturbances.  Echo  signals  arriving 
from  directions  other  than  the  programmed  beam  direction  will 
be  wrongly  identified  in  the  ionograms,  depending  upon  which 
antenna  beam  is  picking  up  these  signals  better.  This  com¬ 
plicates  the  interpretation  of  the  ionograms.  Each  amplitude 
value  in  the  ionograms  is  accompanied  by  a  status  number,  as 
shown  in  the  playback  from  magnetic  tape  in  Figure  1.  The 
insets  in  Figure  1  show  the  optically  weighted  font  (Patenaude 
et  al,  1973)  that  is  used  to  print  the  ionograms.  The  meaning 
of  the  status  values  for  different  time  periods  is  explained 
in  Table  1. 
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An  ionogram  is,  in  essence,  a  snapshot  of  the  iono¬ 
sphere,  and  to  study  the  morphology  of  the  ionosphere  requires 
the  analysis  of  many  ionograms.  The  University  of  Lowell  in 
cooperation  with  AFGL  has  developed  computer  software  to  com¬ 
press  each  ionogram  into  a  linear  array  by  either  integrating 
the  amplitudes  along  the  height  or  the  frequency  axis.  In  tf 
first  case  the  result  is  the  so-called  frequency  characteris¬ 
tic,  in  the  second  the  height  characteristic  (Bibl,  1956; 
Buchau  et  al ,  1978).  For  the  frequency  characteristic  the 
height  intervals  above  and  below  160  km  are  processed  inde¬ 
pendently  giving  the  E  (Figures  2  and  3)  and  F-region  (Fig¬ 
ures  4  and  5)  frequency  characteristics.  The  signal  ampli¬ 
tudes  as  function  of  frequency  for  each  ionogram  are  printed 
as  one  vertical  line.  Because  of  the  photographic  reduction 
the  individual  numbers  become  too  small  to  be  identifiable. 

As  Figure  1  shows,  the  Optifont  (Bibl,  1974)  creates  a  quali¬ 
tative  image,  however.  The  variation  of  the  critical  fre¬ 
quencies  of  the  E  and  F  layers  as  function  of  time  can  con¬ 
veniently  be  read  from  these  characteristics  and,  as  we  will 
show  later,  they  help  in  identifying  periods  of  disturbed 
condit ions . 

Figures  6  and  7  show  height  characteristics  of  the 
E  and  F  layers.  The  dark  traces  (large  integrated  amplitudes) 
indicate  the  approximate  minimum  virtual  heights  of  the  layers 
Since  the  echo  signatures  in  the  vertical  and  backscatter  ion¬ 
ograms  are  quite  different,  separate  characteristics  were  pre¬ 
pared  for  the  backscatter  data  (Figures  8  to  11). 

In  an  effort  to  compare  the  local  magnetic  activity 
with  the  ionogram  data  we  prepared  ionospheric  character ist ics 
for  the  months  January,  April,  July  and  September  1980  (see 
Figures  2  to  41)  covering  four  different  seasons.  Three  char¬ 
acteristic  magnetic  indices  DY,  DZ  and  K  were  derived  from  the 


•  e 


-  i.e 

(*«*) 


Figure  2 


9 


Figure  11 


E  REGION  FREQUENCY  CHARACTERISTICS  »  VERTICAL  IONOGRAMS 
1-16  APR  1980  GOOSE  BAY.  LABRADOR 


k'  -  .  Itf  .  -U  fcii  ■  ■  %  f. 

■'•faP-v  ;  *-i  •  -4  ■+$  ». 


<  <  <  <  « 


.  <  <  •  . 


•  .  ;T 


(•■•1  'l#* 
11 12  AST 


12  AST 


ul 


12  AST 


Figure  12 


Figure  21 


F  REGION 


INTEGRATED  HEIGHTS"  BACKSCATTER  IONOGRAMS 
17-31  JULY  I960  GOOSE  BAY.  LABRADOR 


Figure  31 


38 


E  REGION  FREQUENCY  CHARACTERISTICS  “  VERTICAL  IONOGRAMS 
AUQ  31-SEPT  13  1080  GOOSE  BAY.  LABRADOR 


l/BB  IH.CAB 


F  REGION 

FREQUENCY  CHARACTERISTICS 
•«  8ACKSCATTER  IONOGRaMS 
GOOSE  BAY,  LABRADOR 
SEPT  13-30  19S0  i 


INTEGRATED  HEIGHTS  y 
BACKSCATTER  IONOGRAM3^ 
GOOSE  BAY,  LABRADOR  ? 
BERT  13-30  1930  i 


,  V 

a' 


q 


i 


Goose  Bay  magnetometer  readings  for  the  horizontal  component 
Y  (magnetic  north)  and  the  vertical  component  Z.  The  minimum 
values,  and  Z^,  and  the  maximum  values,  Y^  and  Z^ ,  for  each 
1-1/2  hour  interval  were  available,  measured  in  reference  to 
fixed  values  YQ  and  ZQ  (YQ  =  13,800  y  and  ZQ  =  55,846  y). 

From  these  values  the  monthly  averages  7T  ,  Z. ,  YD  or  Z_  were 
subtracted  for  each  time  period  to  eliminate  the  daily  varia¬ 
tion  and  to  detrend  the  local  K  index  (Mayaud ,  1980). 


Y*  =  Yt  -  7.  ;  Z*  =  Zt  -  I.  ;  Y»  =  Yp  -  Yp ;  Z' 


Zr>  ~  Zn  . 


Magnetic  ranges  were  determined  for  a  three  hour  time  window 
which  was  shifted  along  the  time  axis  in  1-1/2  hour  incre¬ 
ments.  This  gives  an  improved  time  resolution  in  comparison 
with  the  internationally  defined  K  value  which  increments  the 
time  window  by  three  hours  (without  any  overlapping).  The 
ranges  for  a  three  hour  window  were  defined  as 

AY  =  Y  -  Y 

a  IRmax  ILmin* 

and,  respectively,  as 

AZ  ~  ZRmax  “  ZRmin* 

YRmax  the  maximum  of  two  adjacent  Y’  values: 


YRmax(t )  =  MaxlYR(t-l);  YR(t)^ 


Similarly: 


ZRmax ( t )  "  MaxlZR(t-l)’  ZR(t)l 


For  the  minima: 


YLmin (t )  =  MinlYL(t-l);  YL(t) 


ZLmin(t)  =  MinlZL(t-l);  ZL(t)l* 
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The  local  K  coefficient  is  specified  in  Table  2  as  function 
of  the  maximum  of  2AY  and  10AZ.  This  local  K  index  is  not 
the  same  as  the  standard  K  index  because  only  data  from  one 
horizontal  component  (in  the  average  magnetic  field  direction) 
was  available. 

Although  the  K  value  is  widely  accepted  for  sorting 
different  classes  of  data  we  were  able  to  better  relate  varia¬ 
tions  in  the  ionospheric  data  to  magnetic  activity  by  looking 
at  the  magnitude  of  the  horizontal  and  vertical  field  compo¬ 
nents.  While  the  K  parameters  are  equally  sensitive  to  posi¬ 
tive  and  negative  magnetic  bay  events  the  Y  values  permit 
clear  distinction  between  them.  In  one  important  local  event 
(11  Nov  1980)  studied  simultaneously  from  an  aircraft  over 
Greenland  and  at  Goose  Bay,  Labrador,  we  had  a  strong  positive 
magnetic  bay  at  Goose  Bay  and  simultaneously  negative  bays  in 
the  Alaskan  magnetometer  chain  recordings.  The  parameters 
DY  =  Y£  +  Y^  and  DZ  =  Z^  +  Z^  have  been  scaled  for  each  1-1/2 
hour  period;  the  scaled  value  of  DY  and  DZ  covers  the  time 
period  between  t  and  t  -  1-1/2  h,  while  K  includes  variations 
from  the  indicated  time  t  to  (t  -  3)  hours.  Our  plan  is  to 
compare  the  local  parameters  with  the  observed  ionospheric 
phenomena  and  the  planetary  Kp  index.  The  three  local  mag¬ 
netic  indices  are  plotted  in  Figures  42  to  49  for  the  four 
selected  months  of  1980.  To  facilitate  comparison  with  the 
ionospheric  characteristics,  the  same  scale  was  used  for  the 
time  axis.  To  compensate  for  the  offset  of  magnetic  data,  DY 
and  DZ  were  shifted  by  4-1/2  hours,  in  the  translation  from  UT 
to  AST  while  the  K  values  recorded  in  Universal  Time  were 
shifted  by  5-1/4  hours. 

The  F-region  frequency  characteristics  show  a  num¬ 
ber  of  days  with  unexpectedly  low  ionization  at  local  after¬ 
noon.  These  days,  listed  in  Table  3,  are  marked  by  enhanced 
absorption  during  the  day,  indicated  by  the  rugged  fmin  in 
the  E-region  and  F-region  frequency  characteristics,  and/or 
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Max  [2Y,  10Z] 

K 

0  -  14 

0 

15  -  21 

1- 

22  -  28 

1 

29  -  35 

1+ 

36  -  43 

2- 

44  -  53 

2 

54  -  66 

2  + 

67  -  83 

3- 

84  -  104 

3 

105  -  130 

3  + 

131  -  161 

4- 

162  -  198 

4 

199  -  241 

4+ 

242  -  291 

5- 

292  -  348 

5 

349  -  413 

5  + 

414  -  486 

6- 

487  -  568 

6 

569  -  660 

6  + 

661  -  762 

7- 

763  -  875 

7 

876  -  1000 

7  + 

1001  -  1138 

8- 

1139  -  1290 

8 

1291  -  1457 

8  + 

1458  -  1640 

9- 

>_  1641 
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Table  2.  Goose  Bay  Interpolated  Local  K-Values 


Bi  17  Nov  82 


1200  AST 
28 

LOCAL 
17-30  JAN 


1200  AOT 

29 


I  ' 

1200  AOT 


30 


MAGNETIC  INDICES 

1980  GOOSE  BAY,  LABRADOR 


Figure  43 


1100  AST 

104 


1000  AOT 

105 


1000  AOT 

too 


LOCAL  MAGNETIC  INDICES 
1-15  APRIL  1980  GOOSE  BAY,  LABRADOR 


Figure  44 


1200  AST 

119 


1100  AOT 
120 


1200  AOT 
121 


LOCAL  MAGNETIC  INDICES 
16-30  APRIL  1980  GOOSE  BAY,  LABRADOR 


Figure  45 


1200  AST 


1200  AST 


1200  AST 


1 


~T  r  T  '  '  ^  T  1  r  1  1  '  f  '  I  I  l  I  '  I  !  1  "I  I  1  T  I  I  I  I  f  |  |  |  I  I  I  ' 

**T  1200  AST  1200  AST 

*96  197  198 


LOCAL  MAGNETIC  INDICES 
02-16  JULY  1980  GOOSE  BAY,  LABRADOR 

Figure  46 


56 


j Oa.'a.'U 


1200  AST 

257 


1200  AST 

258 


1200  AST 


259 


LOCAL  MAGNETIC  INDICES 
AUG  31 -SEPT  15  1980  GOOSE  BAY,  LABRADOR 

C  MPT  a  LIFT  ovr  ) 


Figure  48 


1100  AST 


the  sharp  decline  of  the  F-region  critical  frequency  and  a 
simultaneous  rise  of  the  virtual  heights.  A  good  illustration 
is  the  comparison  of  a  '’normal”  day,  January  7,  198  0  ,  with 
January  3,  1980,  where  foF2  decreased  rapidly  to  less  than  3 
MHz  at  1700  AST  and  the  F-layer  height  increased  (Figure  50). 
On  January  7,  on  the  other  hand,  the  F  ionization  minimum  is 
reached  at  2019  AST  (Figure  51).  The  relationship  between 
increased  absorption  by  enhanced  particle  dumping  into  the 
lower  ionosphere,  the  formation  of  aurora  E  and  the  vertical 
motion  of  the  F-region  ionization  must  still  be  studied. 

The  new  ionogram  programs  will  certainly  facilitate 
the  description  of  these  important  phenomena  of  the  auroral 
and  subauroral  ionosphere.  On  February  4,  1983,  a  very  dis¬ 
tinct  event  of  the  blowing -out  of  the  F-region  ionization  has 
been  monitored  by  chance  in  Lowell  and  was  visible  even  in 
Digisonde  128  recording  at  Cape  Canaveral,  Florida.  For  one 
hour  the  F-region  height  increased  there  substantially,  while 
h'F2  was  twice  as  large  as  normal,  for  several  hours  in 
Lowell,  Massachusetts.  During  this  time  and  many  hours  there¬ 
after  a  complete  blackout  alternated  with  strong  aurora  E  and 
Es  occurred  in  Goose  Bay,  Labrador.  Although  these  events  are 
much  less  pronounced  during  the  summer  season  we  hope  to  catch 
more  of  the  daytime  F-ionization  breakdown  after  the  recent 
installation  of  a  direct  warning  system,  triggered  by  Goose 
Bay  events.  A  similar  'blow-out'  of  the  F-region  ionization 
will  be  described  in  the  next  section. 

When  the  magnetic  characteristics  are  checked  for 
signatures  one  can  see  that  on  most  of  the  selected  days  the 
horizontal  component  Y  shows  a  marked  increase  starting  be¬ 
tween  13  and  15  h  local  time.  The  onset  of  the  rapid  decrease 
of  the  F-region  ionization  and  the  rise  of  the  layer  as  pre¬ 
sented  in  Table  3  in  AST,  are  indicated  as  arrows  in  the  mag¬ 
netic  charts  of  Figures  42  to  49.  Even  the  average  DY  behav¬ 
ior,  produces  by  overlaying  the  magnetic  data  synchronized 
with  the  ionospheric  event  in  the  superposed  epochs  method , 


shows  a  distinct  positive  bay  of  the  horizontal  magnetic  com¬ 
ponent,  with  the  steepest  slope  at  the  onset  of  the  ionospheric 
event  (Figure  52).  Much  work  remains  to  be  done.  We  plan  to 
investigate  the  chronological  order  of  occurrence  of  the  ir¬ 
regular  magnetic  variation,  the  F  ionization  change  and  the 
formation  of  aurora  type  E  layers. 

The  following  significant  statements,  resulting  from 
the  comparison  between  magnetic  and  ionospheric  data  for  the 
four  months  of  1980  Goose  Bay  data,  can  be  made: 

1.  As  Figure  53  shows,  the  time  of  occurrence  of  positive 
bays  with  average  (Y£  +  Y^)/2  raaxirau111  excursions  larger 
than  25  [y]  is  limited  between  1500  to  0000  Local  Mean 
Time  (AST ) . 

2.  The  time  of  occurrence  of  negative  bays  with  the  same 
maximum  excursions  is  limited  between  2230  and  0730  AST. 

3.  All  significant  positive  bays  are  coinciding  with  fast 
upward  motion  of  the  F  region  ionization  resulting  in 
local  ionization  depletion. 

4.  Most  minima  in  ionization  are  related  to  increases  of 
the  horizontal  field  components. 

5.  Positive  bays  are  related  to  local  precipitation,  while 
the  negative  bays  might  be  related  to  events  occurring 
in  a  latitude  sector  of  90°  to  180°  towards  west. 

Automatic  processing  of  the  Digisonde  ionograms 
facilitates  the  interpretation  in  terms  of  geophysical  phe¬ 
nomena.  A  summary  of  such  techniques  is  given  in  Scientific 
Report  No.  1  'Processing  of  Digital  Ionograms'  by  B.  W. 

Reinisch  (1981). 


AVERAGE  CHANGE  OF  HORIZONTAL  MAGNETIC  FIELD 
COMPONENT  FOR  17  EVENTS  OF  MINIMUM  F  REGION 
IONIZATION  DURING  DAYTIME 


On  11  November  1980,  AFGL’s  Airborne  Ionospheric 
Observatory  (AIO)  onboard  a  KC-135  airplane  monitored  a  day¬ 
time  collapse  of  the  F-region  ionization  over  the  southern 
tip  of  Greenland  (67°  CG  Latitude).  The  aircraft  observations 
were  combined  with  the  ground-based  ionosonde  measurements  at 
Goose  Bay,  Labrador  (65°  CGL),  St.  John’s,  Newfoundland  (58° 
CGL),  Narssarssuaq  (69°  CGL)  and  Godhaven  (76°  CGL)  in  Green¬ 
land.  The  results  of  this  study  were  presented  at  the  XXth 
URSI  General  Assembly  in  Washington,  D.  C.,  August  1981  (K. 
Bibl  et  al,  1981).  Digisonde  ionograms  with  amplitude  (bot¬ 
tom)  and  Doppler  (top)  recorded  by  the  AIO  at  1736  UT  and 
1806  UT  (approximately  1500  LT)  are  shown  in  Figures  54  and 
55.  The  F-layer  critical  frequency  decreased  from  11.2  MHz 
(1.5  x  io12  el/m3)  to  2.7  MHz  (0.9  x  1011  el/m3),  while  the 
minimum  virtual  height  increased  from  235  to  410  km.  During 
this  time,  the  aircraft  was  zigzagging  between  68°  and  73° 

CGL  (Figure  56),  and  changed  local  magnetic  time,  but  returned 
to  almost  the  same  geographic  location.  In  order  to  interpret 
The  aircraft  observations  in  terms  of  the  larger  scale  iono¬ 
spheric  structure  of  the  auroral  dusk  to  midnight  region  we 
used  available  ground-based  ionogram  data  from  stations  lo¬ 
cated  between  65°  CGL  (Goose  Bay)  to  76°  CGL  (Godhaven).  The 
best  information  was  available  from  Goose  Bay  where  12  verti¬ 
cal  and  12  backscatter  ionograms  were  recorded  per  hour  dur¬ 
ing  the  aircraft  mission.  A  quick  check  of  the  hourly  values 
for  foF2  had  not  shown  any  unusual  behavior  for  11  November 
1980.  But  the  F-region  frequency  characteristic  for  this  day 
(Figure  57)  exhibits  a  sharp  decline  in  foF2  starting  at  1510 
AST  (1910  UT ) ,  although  the  decrease  is  somewhat  obscured  by 
oblique  echoes.  The  increase  in  virtual  range  from  about 
270  km  to  400  km  within  15  minutes  is  clearly  visible  in  the 
h*  characteristic  (Figure  58).  The  ionogram  at  1504  AST 
(1904  UT)  shows  three  distinct  traces  clearly  distinguished 
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by  different  Doppler  frequencies  (the  Doppler  ionogram  is  at 
the  bottom  of  Figure  59,  the  amplitudes  are  on  top).  Positive 
Dopplers  indicate  approaching  and  negative  Dopplers  receding 
formations.  The  lowest  trace  in  the  ionogram  has  strongly 
negative  Dopplers,  the  middle  one  slightly  positive  Dopplers 
and  the  highest  trace  very  positive  Dopplers.  This  is  the 
typical  signature  of  a  moving  ionization  hole  the  center  of 
which  has  almost  reached  the  station  (Bibl,  1974). 

The  scaled  peak  plasma  densities  measured  at  the 
different  locations  and  times  where  mapped  in  a  polar  cor¬ 
rected  geomagnetic  coordinate  system  (Whalen,  1970).  Inter¬ 
polation  between  the  discrete  observations  resulted  in  the 
iso-density  contour  plot  shown  in  Figure  60.  The  mid-latitude 
trough  with  plasma  frequencies  of  3  to  4  MHz  extends  from  dusk 
to  midnight.  A  localized  F-region  bite-out  occurs  at  1620 
corrected  geomagnetic  local  time  with  very  steep  gradients  at 
the  western  wall.  A  quiet  reference  day,  26  January  1981, 
shows  no  afternoon  bite-out  (Figure  61). 

If  the  trough  is  believed  to  be  the  result  of  the 
superposition  of  the  corotational  eastward  convection  and  the 
magnetospheric  polar  cap  convection  (Knudsen,  1974;  Spiro 
et  al,  1978),  then  changes  in  the  dawn  to  dusk  magnetospheric 
E-field  will  strongly  affect  the  shape  and  location  of  the 
trough.  However,  we  do  no.  expect  the  stagnation  point  of 
the  convection  to  be  at  local  afternoon  (1600  CGLT)  when  the 
F-region  is  still  illuminated  by  the  sun.  It  is  more  likely 
that  the  hole  is  not  a  direct  part  of  the  premidnight  trough. 
The  observed  increase  in  virtual  height  suggests  an  upward 
motion  of  the  F  ionization  in  a  localized  area,  but  it  is  not 
clear  what  physical  process  would  explain  this  field  aligned 
plasma  flux. 
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2.4  Routine  Ionoeram  Scalin 


The  hourly  vertical  lonograms  of  Goose  Bay  have  been 
scaled  using  URSI  regulations  as  documented  by  Piggott  and 
Rawer  in  the  UAG  Reports  23  (1972),  23A  (1978)  and  50  (1975) 
issued  by  the  World  Data  Center  A,  Boulder,  Colorado.  Since 
these  ionospheric  background  data  serve  as  reference  for  many 
special  studies  we  have  scaled  a  complete  set  of  ionospheric 
parameters . 


Frequencie  s : 


f  oF2 


foFl 


foEs 


Type  Es 
f  oEs2 
Type  Es2 
fbEs 


f  rEs 

MUF (3000)F2 


MUF (3000)Fi 


f  min 

foFOB 

fminOB 


critical  frequency  of  F2  layer 
critical  frequency  of  FI  layer 
critical  frequency  of  E  layer 
critical  frequency  of  Es  layer 

critical  frequency  of  second  Es  layer 

critical  frequency  of  blanketing  Es 

critical  frequency  of  transitory  layer  between 
E  and  F 

cusp  frequency  separating  E  and  Es 

Max.  usable  frequency  for  F2  layer  propagation 
on  3000  km 

Max.  usable  frequency  for  FI  layer  propagation 
on  3000  km 

lowest  frequency  that  shows  echo 
max.  frequency  of  oblique  F  echo 
min.  frequency  of  oblique  F  echo 


Virtual  heights  (ranges): 


h  *  F2 
hpF2 


h’Es 


minimum  height  of  F2  layer 

true  peak  height  if  F  layer  has  parabolic  shape 

minimum  height  of  FI  layer  (if  present)  or  F2 
layer 

minimum  height  of  E  layer 
minimum  height  of  Es  layer 


h’Es2  minimum  height  of  second  Es  layer 

h'FOB  range  of  oblique  F  echo 

The  oblique  sporadic  E  and  F  traces  are  scaled  to  monitor  ir¬ 
regularities  moving  over  the  station.  Hourly  values  are  not 
sufficient  to  analyze  these  irregularities,  but  at  least  they 
indicate  the  disturbed  time  periods. 

Ionograms  from  January  1979  to  September  1982  were 
scaled,  the  values  were  computer  tabulated  and  stored  on  punch 
card  and  on  magnetic  tape.  Actually  two  magnetic  tapes  are 
updated  every  month  by  adding  new  data.  One  tape  is  stored 
at  AFGL,  the  other  at  the  University  of  Lowell.  Data  from  be¬ 
fore  1979  are  currently  archived  only  on  punch  cards,  and 
transfer  onto  magnetic  tape  will  be  done  in  the  future.  Table 
4  shows  the  program  required  to  read  the  tape.  The  computer 
program  that  sorts  and  error  checks  the  scaled  data  also  cal¬ 
culates  the  monthly  medians.  The  median  curves  for  foF2, 
foFl,  foE,  foEs,  h’F2,  h*F,  h*E,  h'Es  and  the  Es  occurrence 
distributions  are  shown  in  Appendix  A  for  the  period  January 
1979  to  September  1982. 

In  winter  1979/80  the  solar  sunspot  cycle  reached 
its  peak.  This  is  evident  by  the  high  foF2  values  between 
1200  and  1500  local  time.  In  the  winter  78/79  the  median 
values  reaches  a  peak  of  12.5  MHz  (February  1979),  in  November 
1979  the  peak  is  14.6  MHz,  in  November  1980  it  is  13.5  MHz  and 
in  December  1981  13.9  MHz.  In  comparison,  the  peak  median 
values  for  the  winter  76/77  was  6.2  MHz  (January  1977)  and  for 
winter  77/78  it  was  9.1  MHz  (February  1978).  The  highest  E 
layer  critical  frequency  occurred  in  June  1980  with  foE  = 

3.9  MHz. 


2.5  Computer  Printing  of  Digisonde  Ionograms 


It  is  often  desired  to  obtain  computer  printouts  of 


83/07/13.  14.17.10. 

PROGRAM  CLAIR 

PROGRAM , CLAIR ( INPUT , OUTPUT , TAPE5= INPUT , TAPE6=0UTPUT , TAPE2 , PUNCH 
1  TAP!  4  , TAPE  1 1 , TAPE1 2 ,  TAPE1 3fTAPE14fTAPE15,TAPE16f  TAPE1 7  ,  TAPE18 , 
1TAPE19,TAPE20fTAPE21 ,TAPE22) 

DIMENSION  IFQF( 12r31 f24>  , IDA( 12 > , X< 31 ) , I AV<24 ) , NCT( 24 ) 

DATA  IDA/31 , 28,31 f 30 , 31 f 30,31 f 31 f 30 r 31 ,30 r 31/ 

NF 1=17 
NF2= 18 
DO  1  K  =  1  , 24 
DO  1  1=1,12 
DO  1  J  =  1 ,31 

1  IFQF ( I , J , K  )  =0 
DO  7  NF=NF1 , NF2 
READ<NF, 18) ID 
WRITE (8, 30) ID 

22  CONTINUE 

READ(NF, 18) ID 
I F ( EOF ( NF ) ) 7 , 8 
8  CONTINUE 

READ ( NF , 1 8 ) I D , I YR , MON , I DA Y , I HR , I FG 
IF ( EOF ( NF ) ) 7 »  32 
32  CONTINUE 

I F ( I DAY . EQ . 0 ) GO  TO  20 

I F ( MON , EG . 0 ) GO  TO  20 

I F  QF ( MON , I  DAY , I HR+ 1 ) =2*IFG 

2  READ(NF, 18) ID 
IF ( EOF ( NF ) ) 7 ,22 

7  CONTINUE 

IM1 =NF1- 10 
IM2=NF2-10 
DO  2  I = I M 1 , IM2 
IM= IDA ( I ) 

IF ( IYR.EQ.7G.AND. I .EG. 2) IM=IM+1 
IF ( I YR . EG . 80 . AND  I .EG. 2) IM=IM+1 
DO  2  J=1 , IM 

WRITE(G,30) IYR, I , J , ( IFQF  < I , J , K ) , K  =  1 , 24 ) 

2  WRITE (4  ) IYR, I , J, ( IFGF( I , J,K ) ,K=1 ,24) 

ENDFILE  4 
REWIND  4 

30  FORMAT ( IX ,301 4) 

18  F  ilRMAT  (  1 2 , 3X  ,  4 1 2 , 43X  ,13) 

19  FORMAT (14) 

STOP 

END 


Table  4 


ionograms  recorded  on  magnetic  tape.  Based  on  existing  soft¬ 
ware  we  developed  a  special  routine  for  this  task. 

The  computer  program  prints  ionograms  from  either 
Digisonde  128  or  Digisonde  128PS  tapes.  For  each  ionogram, 
the  program  prints  a  header  label  with  preface  information,  a 
height  scale  in  km  and  amplitude  data  for  each  frequency.  For 
each  frequency  line,  the  program  prints  the  actual  frequency 
(in  0.1  KHz  steps)  and  126  amplitude  values  (in  4  dB  incre¬ 
ments)  followed  by  the  noise  level  and  the  gain.  This  is  sim¬ 
ilar  to  the  on-line  printing  format  in  the  aircraft.  The  user 
may  specify  whether  the  program  will  print  all  the  ionograms 
on  the  tape,  only  ionograms  at  hourly  intervals  or  only  ver¬ 
tical  ionograms.  The  program  is  written  in  FORTRAN  using 
Control  Data  Corporation’s  FORTRAN  Extended  Version  4  and  will 
need  modification  if  used  on  a  non-CDC  computer.  Figure  62 
shows  an  example  of  such  printout,  photographically  reduced 
to  report  size.  The  print  program  has  been  implemented  on 
the  computer  of  World  Data  Center  A,  Boulder,  Colorado,  where 
it  is  available  to  provide  on  request  (Conkright,  1983)  print¬ 
outs  of  digital  ionograms  from  Goose  Bay.  The  total  magnetic 
tape  data  base  (Digisonde  128  and  128PS  ionograms  since  January 
1975)  has  been  transferred  to  WDC  A. 

A  full  description  of  the  program  is  given  in  an 
AFGL  Research  Paper  (Shirley  et  al,  1983). 


3.0  GROUNDBASED  AND  AIRBORNE  IONOSONDE  OBSERVATIONS 


IN  SUPPORT  OF  THE  OTH  RADAR  TESTS  1980/81 

The  AFGL  airborne  ionospheric  observatory  was  carry¬ 
ing  out  a  number  of  missions  in  1980/81  serving  as  target  for 
the  OTH  radar  during  the  testing  of  the  Experimental  Radar 
System  (ERS)  and  at  the  same  time  collecting  ionospheric  data. 
Vertical,  bistatic  (propagation)  and  backscatter  ionograms 
were  made  at  Goose  Bay,  the  aircraft  and  at  the  OTH  site  in 
Maine.  The  aircraft  has  only  one  HF  antenna  which  had  to  be 
multiplexed  between  the  different  tasks.  The  sounding  sched¬ 
ule  that  was  implemented  is  shown  in  Figure  61,  listing  the 
transmission  and  receiving  times  for  the  Goose  Bay  Digisonde, 
the  aircraft  Digisonde,  the  aircraft  Barry  sounder,  and  the 
aircraft  HF  transponder.  The  aircraft  Digisonde  made  six  ver¬ 
tical  ionograms  per  hour,  and  the  radiation  was  received  as  a 
propagation  ionogram  at  Goose  Bay  operating  its  Digisonde  in 
the  receive-only  mode.  The  Goose  Bay  Digisonde  transmitted  16 
ionograms  per  hour,  ten  of  which  were  received  at  the  air¬ 
craft.  The  Goose  Bay  B  ionograms  used  vertical  transmit  and 
receive  antennas  to  make  vertical  ionograms.  In  parallel,  a 
log-periodic  rotatable  antenna  was  activated  to  transmit  in 
the  direction  of  the  aircraft.  The  C  ionogram  was  programmed 
for  backscatter  observations  using  transmit  and  receive  an¬ 
tennas  pointing  toward  magnetic  north. 

For  the  radar  operators  at  Maine  up-to-date  informa¬ 
tion  on  the  ionosphere  in  the  ray  path  reflection  region  is 
important  for  frequency  management  and  for  radar  range  to 
ground  range  conversion.  To  feed  complete  ionograms  through 
a  dedicated  telephone  line  to  Maine  ULCAR  had  built  the  Iono¬ 
gram  Communicator,  ICOM  (Reinisch  and  Bibl,  1981).  Different 
ionograms  can  be  selected  for  telephone  transfer  by  operating 
ICOM  in  Mode  12,  13,  14  or  15  (Figure  63),  allowing  the  radar 
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Figure  63.  HF  Operations  During  OTH  Missions 
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operator  to  consult  vertical,  backscatter  and  propagation 
(aircraft  to  Goose  Bay)  ionograms. 

Dr.  R.  Penndorf  analyzed  a  large  amount  of  the  prop 
agation  data  and  prepared  a  number  of  internal  reports  sub¬ 
mitted  to  the  Contract  Monitor.  An  AFGL  in-house  test  report 
summarized  the  results:  "AFGL  Airborne  Ionospheric  Observa¬ 
tory,  ERS  Test  Report  No.  2,"  Technical  Memorandum  45,  26 
January  1981  by  J.  Buchau,  E.  Weber,  R.  Gowell,  P.  Krukonis, 
J.  Moore,  K.  Bibl  and  B.  Reinisch. 


4.0  IONOSPHERIC  DRIFT  OBSERVATIONS  AT  GOOSE  BAY 


Motions  of  the  ionosphere  can  be  detected  and  ana¬ 
lyzed  by  means  of  the  Doppler  shift  imposed  on  the  radio  wave 
by  the  moving  reflector  (Pfister  and  Bibl,  1972).  The  Digi- 
sonde  128PS  (Bibl  and  Reinisch,  1978)  in  Goose  Bay,  when  oper 
ating  in  the  Doppler-Drift  mode,  measures  the  full  Doppler 
spectrum  of  the  echo  signals  received  at  each  of  the  four  re¬ 
ceiving  antennas.  Cross-correlating  the  complex  spectra  (am¬ 
plitude  and  phase)  of  the  antenna  signals  results  in  sky  maps 
(Bibl  et  al,  1975)  which  give  the  locations  of  the  reflection 
points  and  the  Doppler  frequency  imposed  by  the  motion  of  the 
reflection  points.  Figure  64  shows  on  the  left  side  the  loca 
tions  of  the  reflection  points  for  four  consecutive  18  second 
measurements  (cases)  labelled  0,  1,  2,  and  3.  The  associated 
Doppler  frequencies  are  displayed  in  the  accompanying  map  on 
the  right.  It  is  important  to  stress  the  existence  of  sev¬ 
eral  reflection  points  for  each  observation.  The  signals  re¬ 
ceived  at  each  antenna  are  the  sum  of  the  waves  reflected  at 
the  different  'source*  points.  It  would  be  meaningless  to 
measure  the  phase  of  the  composite  signal  since  frequently 
the  different  sources  have  comparable  amplitudes.  Simple 
phase  triangulation  techniques  for  the  determination  of  the 
source  location  work  well  only  if  only  one  source  exists. 
Since,  in  general,  the  ionosphere  is  moving  each  reflection 
point  has  its  characteristic  Doppler  frequency: 

d  =  r*  k*v 

TT 

where  k  is  the  wave  vector  pointing  from  the  source  to  the 
receiver  and  v  the  velocity  vector  of  the  reflector.  The 
spectral  analysis  thus  separates  the  signal  contributions 
from  the  different  sources,  and  phase  comparisons  between  the 
antenna  signals  can  be  made  for  each  spectral  component. 


If  a  uniform  drift  velocity  is  assumed  one  can  find 
the  velocity  vector  v  that  best  reproduces  the  observed  Dop¬ 
pler  frequencies  for  a  given  sky  map.  We  have  developed  the 
computer  software  to  the  point  that  large  data  tapes,  with 
drift  data  recorded  at  Goose  Bay,  can  be  processed  efficiently 
outputting  the  three  dimensional  drift  velocity  vectors  as 
function  of  time.  Figure  65  is  an  example  of  the  results  ob¬ 
tained  for  an  eleven  hour  observation  period  from  1800  AST  on 
26  January  to  0500  AST  the  next  morning.  The  azimuth  (bottom) 
and  the  magnitude  (top)  of  the  horizontal  component  of  the 
drift  velocity  are  shown  as  solid  lines.  Each  drift  observa¬ 
tion  was  conducted  simultaneously  at  three  frequencies  and  at 
the  virtual  heights  of  the  corresponding  plasma  densities. 

The  velocities  at  all  three  ranges  are  represented  by  numbers 
and  letters  that  indicate  the  ranges  according  to  the  equation 

r  =  200  +  10X,  X  =  1,  2,  3,  ...,  9,  A,  B,  ... 

For  example  between  2147  and  2200  AST  the  range  indicators 
are  F,  G  and  J  corresponding  to  ranges  350  km,  360  km  and  390 
km.  In  general,  the  different  frequencies  observe  similar 
velocities.  The  solid  line  is  the  average  of  the  three 
observations. 

A  very  systematic  behavior  of  the  drift  direction 
is  seen  in  the  azimuth  display.  Before  local  midnight  the 
drift  direction  is  predominantly  westward,  and  after  midnight 
it  is  mainly  eastward.  A  similar  transition  is  observed  for 
January  20/21  (Figure  66).  The  azimuth  variations  are  signi¬ 
ficant  after  midnight,  when  the  drift  speed  occasionally  sur¬ 
passes  200  m/s.  Before  midnight  the  speed  remains  generally 
below  150  m/s  with  one  exception  around  2000  AST  on  26  January. 
The  vertical  velocity  component  is  shown  in  the  upper  panels 
as  +  and  -  symbols  indicating  upward  and  downward  velocities. 
Most  of  the  time  the  vertical  velocity  is  smaller  than  40  m/s, 
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but  values  up  to  100  m/s  are  observed.  The  observed  hori¬ 
zontal  velocities  agree  in  magnitude  and  overall  directions 
with  the  statistical  results  recently  reported  by  Oliver 
et  al  (1983)  from  Millstone  Hill  Incoherent  Scatter 
Measurements. 

The  midnight  drift  reversal  seems  to  suggest  that 
the  Goose  Bay  F  region  is  part  of,  or  at  least  is  coupled  to, 
the  polar  cap  two  cell  convection  pattern  (Spiro  et  al,  1978) 
Drift  measurements  from  January  82,  January  83  and  March  83 
are  currently  analyzed  to  study  the  midnight  reversal.  The 
observational  techniques  and  the  methods  of  analysis  are  de¬ 
scribed  in  Scientific  Report  No.  6  (Dozois,  1983). 


5.0  POLAR  CAP  IONOSPHERIC  STUDIES  AT  THULE 


In  January  1982  and  January  1983  AFGL’s  KC-135  air¬ 
craft  was  parked  at  Thule  AF  Base,  Greenland  (86°  CGL)  to 
conduct  radio  and  optical  observations  of  the  winter  polar 
cap  ionosphere.  The  Digisonde  operated  at  the  aircraft  HF 
antenna.  Since  no  receiving  antenna  array  was  deployed  we 
could  not  make  incidence  angle  measurements  on  echoes  received 
from  F-region  arcs  and  ionization  patches.  But  all-sky  pho¬ 
tometer  images  (Weber  et  al,  1977),  taken  at  2  to  5  minute 
intervals,  identified  the  locations  of  the  irregularities 
that  were  seen  on  the  Digisonde  ionograms.  The  ionograms 
were  recorded  at  1.5  or  2.5  minute  intervals.  The  large  dy¬ 
namic  effects  in  the  winter  polar  cap  F  region  are  illustrated 
by  the  height  characteristic  of  22  January  82  (Figure  67). 

All  160  hours  of  Digisonde  ionograms  recorded  during  the  Jan¬ 
uary  82  mission  were  processed  with  a  velocity  filter,  based 
on  the  measured  Doppler  and  the  operational  frequency,  and 
presented  as  range/time  characteristics;  an  example  is  shown 
in  Figure  68.  For  instance,  at  23:30  UT  an  F  region  irregu¬ 
larity  approaches  the  station  with  a  horizontal  velocity  of 
approximately  250  m*sec-1,  since  the  characteristic  labelled 
250  m»sec-1  contains  the  highest  signal  amplitudes  of  the 
echoes  trace  that  decreases  in  range  as  function  of  time.  An 
inspection  of  the  range  change  time  history  also  reveals,  that 
the  trace  best  seen  around  23:30  UT  in  the  240  msec-1  velocity 
characteristic  corresponds  to  a  horizontal  bulk  motion  of 
-250  msec-1,  i.e.  the  bulk  velocities  and  the  simultaneous 
Doppler  velocities  are  identical.  A  description  of  the  data 
processing  techniques,  a  summary  of  the  observations  and  the 
conclusions  are  contained  in  the  AFGL  Environmental  Research 
Paper,  No.  792,  by  J.  Buchau,  B.  Reinisch,  E.  Weber  and  J. 
Moore,  "Structure  and  Dynamics  of  the  Winter  Polar  Cap  F  Re¬ 
gion."  The  abstract  of  this  paper  is  reproduced  here. 
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ABSTRACT 


All-sky  photometer  images  and  ionospheric  soundings 
taken  at  Thule,  Greenland  (86°  CGL)  in  December  1979  and  Jan¬ 
uary  1982  reveal  three  groups  of  forms  important  in  the  organ¬ 
ization  of  the  winter  polar  cap  ionosphere. 

1)  The  most  prominent  features  are  sun-aligned, 
generally  unstructured,  subvisual,  F-region  arcs,  extending 
for  more  than  1200  km  (limit  of  all-sky  camera  field  of  view) 
across  the  polar  cap.  These  arcs  usually  drift  from  dawn  to 
dusk  at  speeds  between  100  and  250  m/sec,  however,  stagnations 
of  arc  drift  and  drift  reversals  have  been  observed.  The  arcs 
are  produced  by  soft  particle  precipitation. 

2)  During  a  magnetically  disturbed  period  the  arcs 
disappeared  and  large  patches  of  enhanced  F-region  ionization 
drifted  as  speeds  of  250  to  700  m/sec  across  the  field  of  view 
in  the  antisunward  direction.  Although  arcs  are  produced  by 
soft  particle  precipitation,  preliminary  results  from  the 
Dynamics  Explorer  satellite  do  not  indicate  any  localized  soft 
electron  precipitation  into  the  patches. 

3 )  On  a  few  occasions  both  forms  were  observed 
simultaneously.  Between  F-region  sun-aligned  arcs  drifting 
from  dawn  to  dusk,  small  patches  of  ionization  were  observed 
moving  at  much  higher  speeds  in  the  antisunward  direction. 

Both  the  arcs  and  the  patches  appear  as  strong  lo¬ 
calized  irregularities  in  the  ionospheric  soundings.  The 
Doppler  information  provided  by  the  Digisonde  128PS  was  used 
to  relate  backscatter  traces  to  individual  arcs  or  patches 
and  to  track  these  features  over  more  than  1500  km.  Many  ob¬ 
servations  suggest  specular  reflection  from  electron  density 
enhancements  associated  with  the  optical  forms,  rather  than 
scatter  from  field-aligned  irregularities  embedded  in  the 
arcs  or  patches.  High  velocities  occurred  more  often  during 
magnetically  active  periods. 


F-region  arcs  are  bands  of  enhanced  ionization,  im¬ 
bedded  in  a  background  ionosphere  with  a  base  height  h'F  of 
approximately  250  km  and  a  critical  frequency  of  about  4  MHz 
(2  x  10s  el/cm3).  The  virtual  heights  in  the  ionograms  did 
not  change  during  transit  of  the  arcs  through  the  zenith. 

During  the  active  periods,  when  the  ant i sunward -moving  patches 
were  observed,  the  background  ionization  dropped  to  less  than 
3  MHz  (10s  el/cm3)  while  the  base  height  (h'F)  moved  up  to 
heights  above  400  km.  The  strongly  ionized  patches  (foF2  more 
than  8  MHz),  however,  were  observed  to  reach  a  minimum  virtual 
range  of  about  250  km  during  the  zenith  transit,  leading  to 
rapid  h’F  fluctuations  in  the  order  of  200  km  within  minutes. 
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6.0  AUTOMATIC  REAL  TIME  IONOGRAM  SCALER  WITH 
TRUE  HEIGHT  ANALYSIS,  ARTIST 

In  March  1983,  the  ARTIST  was  installed  at  the 
Goose  Bay  Ionospheric  Observatory  and  connected  to  the  Digi- 
sonde  128PS.  Originally  conceived  for  the  Digisonde  256,  the 
ARTIST  has  been  adapted  to  the  128PS  and  is  now  routinely 
scaling  the  Goose  Bay  digital  ionograms  in  real  time.  After 
many  years  of  development  we  finally  svcceeded  in  complete 
automatic  scaling  of  Digisonde  ionograms  in  real  time.  The 
difficult  task  was  to  scale  ionograms  for  disturbed  condi¬ 
tions.  In  Goose  Bay  about  50%  of  the  ionograms  show  spread  F. 
Goose  Bay  is  a  subauroral  station  and  the  ionograms  display  a 
large  variety  of  features:  quiet  and  disturbed  daytime  re¬ 
cordings,  spread  F  during  the  night,  the  mid-latitude  trough 
moving  over  the  station,  fast  variation  of  the  ionospheric 
parameters  and  frequent  absorption  events.  Four  months  of 
ionograms  for  the  high  sunspot  year  1980  served  as  data  base 
for  the  evaluation  of  the  scaling  and  true  height  algorithms. 
Some  8000  ionograms  for  January,  April,  July  and  September 
1980  were  processed.  The  autoscaled  parameters  for  the  2000 
hourly  ionograms  during  these  periods  were  compared  with  the 
manually  scaled  values.  The  results  of  this  comparison  are 
given  in  Scientific  Report  No.  7  (Reinisch  et  al,  1983). 

Table  5  shows  that  for  89.6%  of  all  ionograms  foF2  was  auto¬ 
scaled  within  ±0.5  MHz,  and  for  95.6%  within  ±1.0  MHz.  This 
breakthrough  in  ionospheric  observational  techniques  will  be 
important  for  communication  prediction  as  well  as  ionospheric 
research.  The  scaling  algorithms  and  the  scaling  accuracy 
are  described  in  Scientific  Reports  Nos.  4,  5  and  7  prepared 
by  Reinisch  et  al  (1982),  Gamache  et  al  (1983)  and  Reinisch 
et  al  (1983)  under  this  contract.  Details  the  true  height 
analysis  are  given  in  two  papers  by  Huang  and  Reinisch  (1983), 
and  Reinisch  and  Huang  (1983). 
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Percentage  of  cases  for  foF2,  MUF3000F  and  M3000F  where  manual  and 
automatic  scalings  fall  within  indicated  limits 
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As  an  8086  Central  Processor  Unit  based  computer, 
the  ARTIST  is  a  Digisonde  subunit,  that  is  housed  in  a 
19”  x  9”  x  22"  chassis.  It  operates  either  in  real  time, 
connected  to  a  Digisonde,  or  it  reads  digital  ionograms  from 
magnetic  tape  (9  tracks,  1600  bps).  Installation  and  opera¬ 
tion  of  the  ARTIST  is  simple:  one  cable  connects  the  ARTIST 
to  the  Processor  of  the  Digisonde,  one  cable  goes  to  the  tape 
drive  and  one  RS232  port  drives  the  printer.  A  CRT  terminal 
is  connected  to  a  second  RS232  port.  The  ARTIST  contains  two 
floppy  disk  drives;  one  drive  contains  the  floppy  with  the 
programs,  the  other  can  be  used  to  store  data  and  to  copy 
disks.  After  activating  the  reset  switch  or  after  power-up 
the  program  is  loaded  from  disk  and  the  program  options  can  be 
selected  from  the  keyboard  by  answering  the  questions  dis¬ 
played  on  the  CRT.  Tables  6a,  b  illustrate  the  procedure. 

It  requires  less  than  1  minute  to  process  one  iono- 
gram  including  the  profile  conversion.  Processing  starts  when 
the  ionogram  scan  is  completed.  During  processing  the  next 
ionogram  can  be  inputted,  this  means  that  a  sequence  of  1  min¬ 
ute  ionograms  can  be  handled.  A  typical  ARTIST  printout  is 
shown  in  Table  6c.  The' top  lines  identify  the  ionogram  and 
give  the  sounding  parameters.  The  next  entry  lists  the  iono¬ 
spheric  parameters:  foF2,  foFl,  h'F,  h'F2,  M3000,  fmin, 
foEs,  MUF (3000),  fminF,  f xl ,  fminE,  foE,  h'E,  and  h’Es.  The 
average  range  spread  (in  km)  for  the  F  and  E  traces  is  shown 
under  QF  and  QE,  and  the  frequency  spread  (in  MHz)  under  FF 
and  FE.  The  complete  F  and  E  traces  are  listed  under  the 
heading  AUTOSCALED  TRACES;  the  virtual  height  (in  km)  is  given 
for  each  sounding  frequency  (0.1  MHz  increments).  These  trace 
data  are  the  input  to  the  true  height  routine.  In  the  near 
future,  the  processed  ionogram  data  will  be  transmitted  di¬ 
rectly  to  AWS  and  the  OTH  radar  to  support  the  radar  range  to 
ground  range  conversion  algorithms  of  the  radar. 


ENTER  INPUT.  PLEASE  AESPOM)  WITHIN  SI  SECONDS, 
OTHERWISE  DEFAULT  VALUE  WILL  BE  CHOSEN. 

EACH  ENTRY  MUST  BE  FOLLOWED  BY  A  CARRIAGE  RETURN. 


•TYPE  1  FOR  ONLINE  PROCESSING  (DEFAULT) 
TYPE  S  FOR  TAPE  READ 

ENTER-2 


•WERE  WAS  TIC  STATION  ? 

SELECT  STATION  NUHBER  FROM  THE  F0LL0WIN6  LIST. 

1— NEW  NATAL,  BRAZIL  DATA 
IX— OLD  NATAL,  BRAZIL  DATA 

6— GOOSE  BAY,  LABRADOR,  DATA  TAKEN  ON  OR  AFTER  JUN  14  1900. 
6X— GOOSE  BAY,  LABRADOR,  DATA  TAKEN  BEFORE  JUN  14  1901. 

0— NONE  OF  TIC  ABOVE  STATIONS 

ENTER-6 


•ENTER  SUNSPOT  NUNBER  (DEFAULT**  FOR  JUN  1903) 
ENTER-* 


•TYPE  1  FOR  PROCESSING  ALL  INPUT 
TYPE  S  FOR  PROCESSING  I  (MORANS 


ENTER— 2 


lONOGRANS(DEFAULT) 
BETWEEN  SPECIFIC  TINE 


ENTER  START  TINEi 

DAY - 74 

HOUR - 16 

NINUTE— 59 

ENTER  END  TIC: 

DAY - 75 

HOUR - 1 

NINUTE— 0 


Table  6a.  ARTIST  Initialization 


•SELECT  PLOT  IONOGRAM  OPTIONS.  ENTER  THE  NUMBER  COS  I  NAT  I  ON  IN  TO  LIST 
BELOU  FOR  THE  DESIRES  PLOTS, 

E6.  ENTER  134  08  123456.  ( DEFAULT- 1  > 

1—  ftJTOSCALED  IONOGRAM 

2 —  OPT I FONT  I0N06RPM 

3— COLOR  I0NO6RRM  (ONLY  PLOT  COMPRESSED  DATA) 

4—  DOPPLER  IONOGAAN 

5—  UNDEFINED  PLOT 

6—  BASELINE  AM)  AMPLITUDE 

ENTER- 


•TYPE  1  TO  PRINT  £4  COMPRESSED  BINS  PER  FREQUENCY  (DEFAULT) 
TYPE  2  TO  PRINT  128  BINS  PER  FREQUENCY 

ENTER-2 


•THERE  ARE  5  CHOICES  OF  PRINT  OUT: 

1—  ULCAR  FORMAT 

2 —  AUS  FORMAT,  PARAMETERS  ft  TRACE  HEIGHTS 

3—  AHS  FORMAT,  PARAMETER  OLY 

4—  AMS  FORMAT,  TRACE  HEIGHTS  ONLY 

5— PRINT  ALL  (DEFAULT) 

ENTER- 


•TYPE  1  TO  REMIND  TAPE (DEFAULT) 

TYPE  2  TO  CONTINUE  READ/MRITE  TAPE. 

ENTER- 


CENTER  2  TO  PRINT  INTERMEDIATE  RESULTS,  OTHERUISE  NO  PRINT. 
ENTER- 


*TO  INTERRUPT  ARTIST  EXECUTION: 

ENTER  "CONTROL  C"  TO  RETURN  CONTROL  TO  THE  OPERATING  SYSTEM 
ENTER  "CONTROL  A"  TO  RESTART  ARTIST  PROGRAM 


INITIALIZATIONS  FINISHED,  START  PROCESSING 


Table  6b.  ARTIST  Initialization 


AUTOSCALED  RESULTS  BY  A.R.T.I.S.T.  ULCAR  JUN  83 

ID  YEAR  DAY  H  M  S  B  E  RHTTQNXZKJG  H-STEP(KN)  F-STEP(KHZ> 
6  1983  74  16:59:  9  88-17  4  2  42  1  8  4  1  7  2  3  5.8  188 


FQF2  F0F1  H'F  H»F2  M3888  FAIN  FOES  MUF  FNINF 

6.8  ***  238,  ««  3.84  1.6  2.3  28.7  2.3 


FXI  FNINEFOE  H'E  H'ES  OF  QE  FF  FE 

7.6  1.6  2.3  128.  128.  *#*«**  t«  .3 


AUTOSCALED  TRACES  [KH1: 


2.  hh»  *****  mh  238.  238.  238.  238.  238.  238.  238. 

3.  238.  238.  238.  238.  238.  238.  238.  238.  238.  243. 

4.  243.  248.  248.  248.  253.  253.  258.  258.  263.  263. 

5.  268.  268.  273.  273.  283.  273.  383.  388.  313.  318. 

6.  323.  328.  338.  351  368.  393.  438.  523.  683. 

1.  M*t*  HW  HHt  tHH  *****  WH  120.  125.  125.  138. 

2.  135.  148.  158.  178. 


TRACE  TRANSITION  FM4Z1  VELM/SI  VEL.  AT  TOPF  CN/S] 
F  -22. 

E  67. 

ES  67. 


F 

E 

ES 

I0N06RAH: 


NORMALIZED  AMPLITUDE  AS  AT  REFLECTION  HEIGHT  188KN  IN  EDB3 
TOPF  2.  1  4.  5.  6. 

F  21  8.  42.  42.  47.  41 


Table  6c.  ARTIST  Output 
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The  next  lines  list  the  Doppler  information  contained 
in  the  echo  trace  and  the  vertical  velocity  component  if  it 
can  be  derived  from  the  Doppler  information.  Since  the  velo¬ 
city  calculations  are  not  final  yet  we  will  not  further  dis¬ 
cuss  them  here.  The  last  entry  gives  the  amplitudes  of  the 
echoes  in  dB;  a  free  space  attenuation  correction  to  100  km  is 
applied  using  the  measured  virtual  heights.  Around  each  Mega¬ 
hertz  the  echo  amplitudes  of  three  neighboring  sounding  fre¬ 
quencies  are  averaged.  These  amplitude  values  can  be  used  to 
correlate  geomagnetic  activity  or  other  parameters  of  rele¬ 
vance  to  radio  absorption. 

For  demonstration  and  verification  it  proved  useful 
to  print  the  full  ionogram  together  with  the  identified  trace 
(defined  by  the  leading  edge  of  the  echo  pulses)  on  the  ARTIST 
output  printer.  In  Figure  69,  the  identified  E-trace  is 
marked  by  a  star  and  the  ordinary  F-trace  by  the  #  symbol. 

The  star  symbol  is  also  used  to  indicate  X-cusp  in  the  F2 
layer.  Numbers,  letters  and  symbols  are  used  to  represent 
the  ionogram.  The  amplitudes  for  the  ordinary  echoes  are 
given  in  multiples  of  10  dB  by  the  numerals  1,  2,  ...»  6 
(positive  Doppler)  and  the  letters  A,  B,  ...»  F  (negative 
Doppler).  X  echoes  are  simply  marked  by  an  X,  and  oblique 
echoes  from  the  north  and  south  by  the  symbols  A  and  V .  A 
software  error  in  the  FORTRAN  operating  system  prevented  us 
from  linking  the  true  height  program  to  the  scaling  program. 

But  this  problem  will  be  solved  soon.  We  will  then  be  able 
to  construct  isodensity  contours  as  a  function  of  time  which 
makes  it  possible  to  study  the  time  development  of  the  elec¬ 
tron  density  profiles.  A  sequence  of  profiles  during  sun¬ 
rise  period  is  shown  in  Figure  70. 
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7.0  CHEMICAL  RELEASE  EXPERIMENTS  AT  NATAL 


During  September  1982  AFGL  conducted  ionospheric 
modification  experiments  near  the  geomagnetic  equator  at 
Natal,  Brazil,  in  a  collaborative  international  program  to 
test  the  current  theories  for  generating  plasma  instabilities 
in  the  equatorial  ionosphere.  Rocket  borne  explosive  re¬ 
leases  of  K20,  C02  and  N2  created  plasma  depletions  at  the 
ionospheric  F-region  ledge.  Barium  releases  produced  deple¬ 
tions  at  electrical  polarization  fields.  ULCAR  participated 
in  these  experiments  with  Digisonde  observations  from  the 
ground  and  from  AFGL’s  KC-135  aircraft.  An  overview  of  the 
different  observations  and  the  preliminary  conclusions  drawn 
are  contained  in  the  Proceeding  of  the  3-5  March  1983  BIME/ 
Coloured  Bubbles  Preliminary  Data  Review  Meeting,  issued  by 
R.  Narcisi,  AFGL,  in  April  1983.  ULCAR  scientists  presented 
or  coauthored  a  number  of  papers  during  this  meeting  (Szusz- 
czewicz  et  al,  1983;  Buchau  et  al,  1983;  Reinisch  et  al, 

1983a;  Buchau  et  al,  1983a;  Reinisch  and  Buchau,  1983). 

In  June  1982,  ULCAR* s  Digisonde  was  installed  at  the 
West  Catre  AFB  (5°  55* S,  35°  15*W,  -6.4°  ML)  near  Natal  to 
collect  background  information  on  ionospheric  conditions  at 
Natal.  A  seven  antenna  crossed-loop  array  was  deployed  per¬ 
mitting  the  identification  of  the  polarization  of  the  echo 
returns  and  their  approximate  angle  of  arrival.  Amplitude, 
angle  of  arrival  and  Doppler  information  made  it  possible  to 
identify  and  track  the  effects  of  the  chemical  releases. 

The  day  to  day  variations  of  the  height  of  the  F 
ledge  is  shown  in  Figure  71  which  displays  the  diurnal  varia¬ 
tion  of  h*  for  4  MHz  (2  x  105  cm“* )  from  13  August  to  18 
September  1982.  The  dashed  section  indicates  the  presence  of 
spread  F.  The  ^Kp  for  each  day  is  given  on  the  right  of  each 
diurnal  curve.  The  four  rocket  releases  are  shown  as  arrows 
on  days  251,  256,  260  and  261.  The  electron  density  profiles 


and  their  changes  have  been  calculated  for  the  sunset  periods 
to  determine  the  rise  velocity  of  the  F-region  ledge.  Figures 
72  through  75  illustrate  the  profile  changes  for  the  four 
launch  days.  Similar  plots  were  prepared  for  the  aircraft 
location  from  the  AIO  ionograms.  For  the  BIME  I  launch  the 
rise  velocity  decreases  from  1.4  km/min  1830  LST  to  0.6  km/min 
at  release  time  (Figure  76).  No  spread  F  development  was 
triggered  by  the  HjO/COj  release.  For  BIME  II  the  rise  veloc¬ 
ity  of  the  F-layer  at  time  of  release  was  1.7  km/min  at  Natal 
(Figure  77)  which  according  to  current  theories  should  permit 
development  of  instabilities  and  spread  F.  However,  no  clear 
indication  for  such  development  was  observed.  By  comparing 
the  Natal  and  AIO  profiles  we  found  a  substantial  east-west 
tilt  of  the  F-region  with  a  gradient  in  ionization  of  about 
9%.  The  rise  velocity  of  the  F-layer  over  the  aircraft,  i.e. 
in  the  neighborhood  of  the  rocket  release,  was  also  smaller 
than  at  Natal. 


log  electron  density 
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Figure  76 


8.0  IONOSPHERIC  HEATING  EXPERIMENTS  AT  ARECIBO 


In  September  1981  AFGL's  Airborne  Ionospheric  Ob¬ 
servatory  participated  in  ionospheric  heating  experiments  at 
Arecibo.  The  HF  heater  antenna  is  located  some  20  km  north 
of  the  Arecibo  Observatory  in  Puerto  Rico.  The  antenna  array 
was  driven  by  transmitters  outputting  up  to  400  kW  of  RF 
power.  Frequencies  of  7.3  and  5.1  MHz  were  selected  for  the 
experiments,  that  were  conducted  at  night.  Optical  and  sounder 
observations  were  made  from  the  aircraft  which  was  shuttling 
back  and  forth  on  a  magnetic  meridian  (approximately)  in  order 
to  map  the  heated  volume  with  a  ray  path  to  a  geostationary 
satellite,  which  was  providing  the  signal  for  250  MHz  scintil¬ 
lation  measurements.  No  scintillations  were  observed,  prob¬ 
ably  because  of  the  low  power  level  of  the  heater's  RF  emis¬ 
sion.  The  optical  instruments  also  did  not  show  airglow  sig¬ 
natures  that  could  be  attributed  to  the  heating.  The  foF2 
values  for  the  five  round-robin  flights  are  shown  in  Figure 
78.  It  was  expected  that  the  HF  heating  would  cause  develop¬ 
ment  of  spread  F.  However,  no  large  scale  spread  F  was  seen 
during  these  nights. 

At  certain  aircraft  locations  (see  tick  marks  in 
Figure  78)  the  Digisonde  observed  oblique  echoes  which  seemed 
to  originate  in  the  heated  volume.  This  was  verified  later 
by  playing  back  the  tape  recorded  ionograms  through  a  Doppler 
filter.  The  oblique  echoes  have  the  correct  Doppler  shifts 
expected  for  waves  reflected  in  the  heated  region  and  received 
on  the  moving  aircraft.  The  aircraft  Digisonde  showed  oblique 
echoes  each  time  when  the  plane  defined  by  the  normal  onto  the 
magnetic  field  direction  and  the  location  of  the  aircraft  was 
intersecting  the  heated  volume.  This  suggests  that  the  heat¬ 
ing  generates  field  aligned  irregularities  (F.A.I.)  which  form 
a  reflecting  grid  for  the  radio  waves. 
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Figure  7  9  shows  a  sequence  of  four  ionograms  re¬ 
corded  during  round-robin  #3  when  the  aircraft  was  crossing 
the  critical  location.  Only  the  second  and  third  ionogram 
(lower  left  and  upper  right)  show  the  echoes,  the  first  and 
fourth  ionogram  has  only  vertical  echoes.  The  ionograms  in 
Figure  79  show  the  Doppler  frequencies  rather  than  the  usually 
presented  amplitudes  of  the  echoes.  It  is  interesting  to  note 
that  the  F.A.I.  echoes  have  zero  Doppler  shift  which  is  ex¬ 
pected  since  the  aircraft  was  flying  in  the  magnetic  meridian 
plane,  i.e.  the  aircraft  velocity  vector  was  at  right  angle 
to  the  wave  propagation  vector.  During  round-robin  #5  the 
aircraft  was  flying  normal  to  this  plane  toward  and  away  from 
the  heated  volume,  and  the  F.A.I.  echoes  showed  alternat ingly 
positive  and  negative  Doppler  frequencies. 

On  September  28,  the  aircraft  flew  over  the  heater 
antenna  array  with  the  Digisonde  operating  in  the  fixed  fre¬ 
quency  mode  at  5.1  MHz  (the  heater  frequency)  to  perform  an¬ 
tenna  calibration  measurements.  The  results  from  these  mea¬ 
surements  are  discussed  in  Scientific  Report  No.  3  (Reinisch 
et  al,  1982). 
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9.0  DIGISONDE  128  INSTALLATION  AT  KENNEDY  SPACE  CENTER 


WITH  REMOTE  DATA  LINK 

A  Digisonde  128  which  had  been  operated  by  AFGL  near 
Keflavik,  Iceland,  had  been  sent  to  Lowell  for  a  thorough 
check-up  and  for  installation  of  a  new  Olivetti  Thermal  Dot 
Printer.  This  printer  is  driven  by  a  microcomputer  built 
into  the  Digisonde  128,  similar  to  the  design  of  the  output 
interfaces  of  the  new  Digisonde  256. 

In  September  1982  the  modified  Digisonde  128  was 
installed  at  Kennedy  Space  Center,  Florida,  replacing  the  old 
prototype  Digisonde  128,  operated  at  this  location  for  sev¬ 
eral  years.  Although  the  Olivetti  printer  did  not  fulfill 
the  expectation  in  reliability,  as  it  had  achieved  in  Belgium 
with  the  new  Digisonde  256,  the  replacement  system  was  a  sub¬ 
stantial  improvement  compared  to  the  old  sounder.  A  tele¬ 
phone  link  for  data  transfer  was  installed  between  Kennedy 
Space  Center  and  the  Air  Weather  Station  at  Patrick  AFB. 

A  remote  printer  provides  all  pertinent  information  of  the 
ionogram  and  facilitates  with  its  frequency  and  range  grid 
and  a  MUF3000  overlay  the  fast  scaling  of  the  most  important 
ionospheric  parameters  by  AWS  personnel.  The  Digisonde  sys¬ 
tem  is  operated  by  RCA  maintenance  personnel  who  also  have 
erected  the  small  vertical  rhombus  transmitter  antenna  and  a 
single  turnstile  loop  antenna  for  reception  following  the 
advice  of  ULCAR.  After  installation  of  the  system  an  intro¬ 
ductory  lecture  on  data  analysis  was  given  for  the  AWS 
personnel. 
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